Effects of electron-electron interaction on the conductance of open quantum dots 
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We study the effect of electron-electron interaction on the conductance of open quantum dots. We 
find that Coulomb interactions (i) do not affect the ensemble averaged conductance (G) if time- 
reversal symmetry has been broken by a magnetic field, (ii) enhance weak localization and weak 
anti-localization corrections to (G) in the absence of a magnetic field, (iii) increase conductance 
fluctuations, and (iv) enhance the effect of short trajectories on the conductivity of quantum dot. 
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The phenomenon of Coulomb blockade in quantum 
dots is commonly associated with dots that are coupled 
to the outside world via tunneling contacts 0. In these 
systems, the total charge on the dot is quantized in units 
of the electron charge e. Once the temperature T be- 
comes smaller than the charging energy Ec (the energy 
needed to add an extra electron to the system), trans- 
port through the quantum dot is suppressed, unless the 
system is tuned to a charge degeneracy point. On the 
other hand, if the connection to the outside world is via 
ballistic point contacts with a conductance much larger 
than the conductance quantum 2e 2 /h, charge is no longer 
quantized, and the Coulomb blockade is lifted. 

In recent years, it has become possible to study the in- 
termediate regime, of quantum dots with ballistic point 
contacts with only a few propagating channels at the 
Fermi level, so that their conductance is ~ e 2 /h. For 
temperatures T comparable to the spacing A of single- 
particle levels in the quantum dot, the conductance of 
these open quantum dots shows mesoscopic fluctuations 
of the same order as the average ||, Theoretically, 
the large mesoscopic fluctuations are understood within 
the framework of random-matrix theory [Q-Q , which em- 
ploys a picture of non-interacting electrons, and does not 
account for Coulomb blockade effects. A partial justifi- 
cation for this approach can be obtained from the work 
of Furusaki and Matveev || , who calculated the effect of 
electron-electron interactions on the conductance G of a 
quantum dot with ballistic single-channel point contacts 
in the limit Ec >T> A, where there are no mesoscopic 
fluctuations. They found G — e 2 /h, indicating that there 
is no Coulomb blockade effect on the conductance for bal- 
listic point contacts in this limit. At the same time, the 
presence of weakly reflecting tunnel barrier in the contact 
will drive the system into a state of Coulomb blockade, 
in which transport is inhibited. 

That this is not the complete picture was explained 
in a recent work by one of the authors and Glazman |jj , 
who showed that as a result of the interplay of mesoscopic 
fluctuations and electron-electron interactions, Coulomb 
blockade effects on the capacitance and the tunneling 
density of states of a quantum dot persist, even at a per- 



fect transparency of the point contacts. This so-called 
"mesoscopic charge quantization" urges us to reconsider 
the effect of Coulomb interactions on the conductance of 
a quantum dot with ideal point contacts. In this letter 
we report a calculation the conductance of such a quan- 
tum dot and its mesoscopic fluctuations, in the presence 
of electron-electron interactions. 

We consider a quantum dot connected to electron 
reservoirs via two leads (labeled 1 and 2), see Fig. [I]. 
In each of the leads there are N propagating channels at 
the Fermi level Ep. For the Hamiltonian of the system 
we take 



Ti. = TLi 



Tic, 



(la) 



where TLf is the Hamiltonian of the non-interacting elec- 
trons and Tic is the interaction Hamiltonian, correspond- 
ing to a capacitive interaction inside the quantum dot, 



Hc = ^(Q-ne)\ Q = e V f drfy a . (lc) 

6 tT=±^ X>Q 




x=0 



FIG. |l|. Quantum dot (dotted) connected to two ideal 
leads (labeled 1 and 2). 

The index a denotes spin. The potential U(f) describes 
the confinement by the gates surrounding the quantum 
dot and the leads. The leads are located at x < 0, and the 
dot at x > 0, see Fig. |l|. As is explained in Ref. M, the in- 
teraction Hamiltonian Tic with the capacitive interaction 
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provides a sufficient description of the electron-electron 
interactions in a metal of semiconductor quantum dot if 
the electrons explore the dot ergodically before they exit 
the dot through one of the point contacts, i.e. if all the in- 
volved energy scales are much smaller than the Thouless 
energy of the dot Et- It is in this regime that the con- 
ductance statistics become universal and random-matrix 
theory can be used as a model for the non-interacting 
Hamiltonian "Hi? 0. 

To summarize the result of our calculation, which is 
described below, we find that the presence of the inter- 
action Hamiltonian (Jlc[) enhances mesoscopic (quantum- 
interference) corrections to the conductance. In particu- 
lar: 

(i) . There is no interaction correction for the en- 
semble averaged conductance (G) in the presence of a 
time-reversal symmetry breaking field (unitary ensem- 
ble). This is different from disordered bulk systems, 
where interactions do have an effect on the conductance, 
even in the presence of a weak magnetic field |Io| . 

(ii) . In the absence of a magnetic field (orthogonal 
ensemble), the weak-localization correction to (G) is en- 
hanced by the interactions. We have computed (G) in a 
perturbation series in NA/T, where N is the number of 
channels in each of the point contacts and A is the mean 
level spacing in the quantum dot, 



N 



c N NA 
~N 8tt 2 T' 



(2) 



Here c/v is a numerical constant ranging from c\ rj 3.18 
to Coo = 7r 2 /6. The conductance is measured in units of 
2e 2 /h. The second term in Eq. (|||) is the usual weak- 
localization correction for non-interacting particles 0, 
the third term is the interaction correction. Within the 
model (1) there is no dephasing, which would give rise 
to a suppression of weak localization. The origin of the 
absence of dephasing is that all interference occurs inside 
the dot, and the interaction Hamiltonian does not pro- 
duce excitations during the course of such motion. In the 
presence of spin-orbit scattering (symplectic ensemble), 
the corrections to the conductance [second and third term 
on the r.h.s. of Eq. (||)] are multiplied by —1/2. Hence, 
in this case, Coulomb interactions enhance the conduc- 
tance. We expect that the conductance distribution satu- 
rates when irT becomes comparable to the level broaden- 
ing NA/-K. Hence for low temperatures, the interaction 
correction to (G) remains small as 1/N, and there is no 
transition to the regime of real Coulomb blockade. 

(iii). The capacitive interaction in the Hamiltonian (1) 
enhances the mesoscopic conductance fluctuations. In 
the unitary ensemble, the conductance G shows sample- 
to-sample fluctuations with variance 

with cat ~ 6.49 for N 3> 1. In the absence of a magnetic 



field, the fluctuations are larger by a factor 2, while the ef- 
fect of strong spin-orbit scattering is to reduce var G by a 
factor 4. The first term on the r.h.s. of Eq. (^), which rep- 
resents the conductance fluctuations for non-interacting 
electrons, is derived in the limit of large channel numbers, 
N 3> 1, and agrees with the result obtained by Efetov 
. The second term is the interaction correction. 

(iv). The above results are for ideal ballistic point 
contacts and for ergodic quantum dots. The presence of 
a weakly scattering impurity in the point contact with 
reflection amplitude r gives rise to a correction to the 
average conductance that behaves as 



N 2 
2tt' 



where c\ « 5.32 and cm — ► 4 for N ^> 1. For N = 1 this 
result has been obtained previously in Ref. [||. However, 
not only an impurity in the contact, but also any other 
direct scattering process that acts on a time scale shorter 
than the Coulomb time t c oc 1/Ec (like scattering from 
short trajectories connecting the point contacts), has the 
same effect on the conductance. The effect of all direct 
processes is summarized by the replacement of the re- 
flection probability \r\ 2 by the difference \r\ 2 — \t\ 2 of the 
probabilities for direct reflection and transmission, re- 
spectively. In particular, we expect that as a result of 
interactions, the conductance of an open quantum dot 
with a direct path connecting source and drain contacts 
increases as the temperature T is lowered. 

Let us now proceed with the formulation of the theo- 
retical framework used for the derivation of these results. 
The conductance is computed from the current-current 
correlator in imaginary time, 



G = — [ dtn(it + — 

H(r) = (T T /(r)/(0)). 



(5a) 
(5b) 



We linearize the spectrum in the leads. In each lead 
there are N channels, labeled j = 1, . . . , N for lead 1 and 
j = N + 1, . . . , 2N for lead 2. The current / is defined in 
the leads, just outside the quantum dot, 



2/V 



3=1 <x=± 



{j J {ja( X )'4 } Lja( X ) ~ ^Rjai^Rjai^ 



(6) 



where Vxja {trjit) anc ^ ^Lja (4>Rja) denote creation and 
annihilation operators for left (right) moving particles 
with spin a in channel j, and v(j) = 1 if 1 < j < N 
and —1 else. (We have chosen units such that the Fermi 
velocity Vp = 1/27T.) We take the limit x f at the end 
of the calculation. 
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The dynamics of the quantum dot is dealt with in the 
effective action method of Ref. || . Following Ref . || , the 
charge Q of the quantum dot in the interaction Hamil- 
tonian Tic is replaced by the charge Ql of the leads. 
Such a replacement is allowed, because the total number 
of particles (in the dot and in the leads) is conserved. 
Next, the degrees of freedom of the quantum dot are in- 
tegrated out, in favor of a "mirror" copy of the leads at 
x > 0, and an effective action S e g, to be defined below. 
The Hamiltonian of the system now reads 



. 2N 

3=1 <7=± 



dx 



nef 



Qi 



2N 

3=1 o-=± 



dx : fl>1 itr if> Lia + V'LV' 



Rj H>Rja 



(7a) 
(7b) 



The scattering of particles from the quantum dot is repre- 
sented by the imaginary-time effective action S e $ acting 
at x = 0, 



<5 e fi - 

X Li4 



1/T 



dn dr 2 0) + ipm a (n; 0)] 



in 



r-i) [ipLj<T< (t 2 ;0) + ipRja' (r 2 ; 0)], 

where ip(r;x) — e HT ip(0',x)e~' HT and Tp(r;x) = 
if>'{—T',x). The kernel Li^ aa i(r) is a hermitian matrix, 
related to the scattering matrix $y.<, a < of the quantum 
dot @, 

L(u n ) = f ' e l ^ T L{ T )dT 



J 



i l-sw i 



A-rti 1 + S(uj n ) 



— sgn(w„), 



(8) 



wherea->„ = (2n+\)irT is the Matsubara frequency. Note 
that in Matsubara representation, the kernel L(u> n ) sat- 
isfies L(u> n ) = L(—U) n y. The current correlator II(t) is 
now calculated as a thermal average with respect to the 
Hamiltonian TL and with the effective action S c $ , 



H(r) = 



(T T /(r)/(0)e' s ^)w 



eff \ 



(9) 



To compute the thermal average of the interacting sys- 
tem defined by Eqs. (@)-@, we first bosonize the one- 
dimensional Hamiltonian TL of Eq. (Q) . In this way the in- 
teraction becomes quadratic in terms of the boson fields, 
and can be dealt with exactly. To account for the scatter- 
ing from the quantum dot, which is represented by the ef- 
fective action S c g , we perform an expansion up to second 
order in powers of the scattering matrix S. To be pre- 
cise, we first expand n(r) in powers of the action <S e ff, and 
then we expand S c s in powers of the scattering matrix S 
using that up to order S 2 we have 2iriL(bj n ) = —S(uj n ) if 
uj n > and S^(— u> n ) otherwise, see Eq. (||). Note that in 
Ref. H the kernel L was used as the expansion param- 
eter. The expansion in the scattering matrix S instead 
of the kernel L guarantees that in the absence of inter- 
actions mesoscopic fluctuations are fully accounted for. 
Only corrections to the conductance that depend on the 
interplay of interactions and mesoscopic fluctuations are 
taken into account perturbatively. The perturbation ex- 
pansion is valid in the regime T 3> NA. We find in this 
way 



G 



Af ttT 1 f°° f 

- - — / dt, dt 2 st kw Ah)Si hw Ah>{mi) - 

Jo k,l,tr,tr> I 



Tsin- 



ds- 



2s + t 2 +h 



sinh[(< 2 - h)-nT\ ' AN 
sinh[(i 2 + ti + s + QttT] sinh[(s - t c )nT] x 1/4W ~ 



sinh[(s + t 2 )irT] sinh[(s + t^irT] \ sinh[(ti + t c )i:T) sinh[(< 2 + t c )irT] 



(10) 



Here t c = ir/2NEc"f, 7 being the Euler constant. The real-time or Lehmann representation of the scattering matrix 
is defined through S(w n ) = / °° S'(t)e~"" t , uj n > 0. 

Equation ( |ici|) contains the conductance of a specific sample in the presence of the interaction (J^) , up to second order 
in the scattering matrix S(t) of the quantum dot. The first term between brackets is nothing else than the Landauer 
formula for the conductance of the non-interacting system. The remaining term is the interaction correction. The 
perturbation theory in powers of the scattering matrix is arranged in such a way, that we handle the non-interacting 
"pole" contributions exactly, and do a perturbation theory in the interacting "branch-cut" terms. The sample-specific 
conductance ( |lfj| ) was obtained by a thermal average of the interacting system. Next, we need to take a mesoscopic 
average over an ensemble of quantum dots. Such an ensemble can be obtained by varying e.g. the shape, impurity 
configuration, or the Fermi energy of the quantum dot. The correlators of the scattering matrix for such a mesoscopic 
ensemble have been studied extensively in the literature [7J. Although the two-point correlator of scattering matrix 
elements is known exactly for our purposes, within the regime T 3> NA where the expansion ( [To| ) is valid, it 
is sufficient to use the asymptotic large- N formulas for these correlators. In the absence of spin-orbit scattering, the 
scattering matrix is diagonal in spin space, S a b- t a<j' = S a b5 aa > . For the unitary ensemble (presence of a time-reversal 
symmetry breaking magnetic field), we have 
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(S: b (t)S a , b ,(t>)) = 6 aa ,5 bb ,-e- NAt ^S(t> t), 
while in the presence of time-reversal symmetry (orthogonal ensemble) we find 



(11a) 



(S: b (t)S a , b ,(t>)) = (5 aa ,5 w + 6 ab ,5 ba ,)-e-^ N+ V At ^6(t> t). 

Z7T 

With spin orbit scattering in the presence of time-reversal symmetry (symplectic ensemble), we have 
(S: b .^ T (t)S a , b ,, a , T ,(t')) = (<W<W<W<W - S&foM^ivyM^e-^-W^Stf - t), 



(lib) 



(11c) 



where a y is the Pauli matrix. For interaction corrections to leading order in NA/T, higher order correlators of S can 
be taken Gaussian. The ensemble average of Eq. ( 10|) with the help of the correlators (|lTJ) yields the weak- localization 
correction (Q) and the conductance fluctuations (|3). (The observation that (G) = N/2 in the unitary ensemble follows 
from symmetry properties of the ensemble of scattering matrices, and holds to arbitrary order in perturbation theory.) 
The theory of Ref. (§], where the effect of weakly reflecting impurities in the contacts was considered, is recovered by 
setting Sij-aa>(t) = rj5(t)8ij5ac!' , r j being the reflection amplitude in channel j. 

For spinless particles (which might e.g. be realized in spin-polarized quantum dots), Eq. (^) remains valid, provided 
one substitutes N — + N/2. However, if the leads contain only one channel (N — 1), an extra rt-dependent term 5G, 
has to be added, 



6G = 



ttT 2 



ds 



dtidt2 



Ree" 



l (2s + ti + t 2 ) [Si 2 (ti)S 2 i(t 2 ) " 5n(ii)5 22 (t 2 )] 



v /sinh[(i c + t^irT] sinh[(i c + t 2 )?rT] sinh[(£ c + s + t 2 + t^irT] sinh[(s - t c )nT] 



(12) 



In this case the weak localization correction to the con- 
ductance is enhanced, 



A t c nT 
In , 

2c ' 



16ttT 



(13) 



where c w 0.36. Eq. ( p_2[ ) bears explicit reference to the 
particle number n, and induces an explicit n-dependence 
in the conductance fluctuations, 



(G(n)G(n')) = 



cA 2 



■cos[27r(ri-ri')] m2 (ic7rT) 



128tt 2 T 2 
+ n-independent terms, 



(14) 



where c = 5/4 (1/2) in the presence (absence) of time- 
reversal symmetry. In our perturbation theory, such an 
explicit n-dependence, which reflects the discreteness of 
charge, is absent for particles with spin. The origin for 
the appearance of this extra term for spinless particles 
and N = 1 is the same as in Ref. ||. For spin 1/2 
electrons, effects of the discreteness of charge appear in 
general in 47Vth order in perturbation theory in S. 

We close with a remark on the issue of dephasing, 
which plays an important role in the experiments on open 
quantum dots Q. While the experiments indicate that 
the dephasing time diverges as T~ l as the temperature 
T — > 0, the theoretical prediction for a closed quantum 
dot is oc T~ 2 |Q. As our model (1), which takes 
the effect of open contacts into account, does not give 
rise to significant dephasing effects, we believe that the 
fact that a quantum dot is open by itself does not lead 
to additional dephasing, and cannot explain the puzzling 
discrepancy between experiment and theory. 
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